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Abstract 

We measure the diffusion rate of Chern-Simons number in the (l+l)-dimensional 
Abelian Higgs model interacting with a realistic heat bath for temperatures be- 
tween 1/13 and 1/3 times the sphaleron energy. It is found that the measured 
rate is close to that predicted by one-loop calculation at the lower end of the 
temperature range considered but falls at least an order of magnitude short of 
one-loop estimate at the upper end of that range. We show numerically that the 
sphaleron approximation breaks down as soon as the gauge-invariant two-point 
function yields correlation length close to the sphaleron size. 



* Supported by the Swiss National Science Foundation. 



1 



Anomalous electroweak baryon-number violation may have played an important role 
in setting the baryon number of the Universe to its present value |l|, U . At temperatures 
above the gauge-boson mass scale electroweak baryon-number nonconservation is domi- 
nated by hopping over the finite-energy barriers separating topologically distinct vacua 
of the bosonic sector. Determination of the corresponding transition rate is a challenging 
nonperturbative problem, even in the range of validity of the classical approximation. At 
the lower temperature end of that range the barrier crossings are likely to occur in the 
vicinity of the saddle point (known as a sphaleron) of the energy functional. The rate 
can then be estimated using a field-theoretic extension of transition-state theory (TST) 
fl7| , |18|j . At higher temperatures this analytic tool is no longer available, and direct mea- 
surement of the rate in real-time numerical simulations of a lattice gauge-Higgs system 
is the only remaining possibility. 

Analytical saddle-point estimates of the rate in the Standard Model are complicated 
by the fact that the corresponding sphaleron field configuration is not known exactly. 
At the same time, numerical real-time simulations of that system in its low-temperature 
regime carry an enormous computational cost and are yet to be performed fl . I n this 
situation lower-dimensional models become a very useful test ground on which activation- 
theory predictions can be confronted by numerical experiments ||, ||, ||, 0, ||, ||, [10|, [11], HI) 
|T5| , |I"5|j . For this reason the (l+l)-dimensional Abelian Higgs model (AHM) studied 
numerically in this work has attracted much attention recently || [|, || [TT[ [12). 

Determination of the transition rate should include its proper averaging over the 
canonical ensemble in the phase space of a system in question. One way to achieve that 
would be to generate the canonical ensemble of initial configurations, subject each of 
these configurations to the Hamitonian evolution, and average the transition rate over 
the initial states of the system. Such procedure, while being perfectly valid, is very costly 
computationally. To date, a single or a small number of initial configurations have been 



used in Hamiltonian simulations of AHM M, 12]. In addition, preparing initial configura- 



tions in case of a gauge theory presents a technical difficulty: if a standard importance- 
sampling method is used, resulting configurations will in general violate Gauss' law. A 
special cooling procedure is required to eliminate static charge ||. It is not clear that 
such cooling does not cause the sample to deviate from the intended canonical ensemble. 

Another way to obtain the canonical ensemble average of the rate is to replace Hamil- 
tonian evolution by evolution in a heat bath. The simplest form of the latter is imple- 
mented using phenomenological Langevin equations of motion. While Langevin approach 
guarantees thermalization of the system, it does so at the expense of introducing an arti- 
ficial viscosity parameter, thereby altering bulk dynamical properties of a field-theoretic 
system. Numerical studies performed on different models show that transition rates in- 
deed strongly depend on viscosity |7|, . Recent analytical work [[H]] has also shown the 
impact of heat-bath properties on quantities of transition-rate type. In case of AHM, 
there also is a technical difficulty with the conventional Langevin approach: in order 
to maintain gauge invariance, one is forced to use polar coordinates for the Higgs field; 
whenever the latter vanishes, the equations of motion are singular 0, [12]. 

Recently we have proposed and tested a new method in which a field-theoretic system 
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interacts with a heat bath at its boundaries ||. The heat bath is constructed so as to 
imitate an infinite extension of the system beyond the boundaries. Technically this means 
that the fields in the bulk of the system evolve according to the Hamiltonian equations of 
motion, while boundary fields are subject to Langevin evolution with a non-Markovian 
friction kernel and colored noise. Our construction approximates a natural situation 
in which open systems are immersed in a similar environment. In this way, dynamical 
evolution and canonical ensemble averaging occur at the same time while bulk dynamical 
properties of the system are intact. Moreover, the new procedure does not suffer from 
the technical difficulties of the two old ones. In this work we apply the realistic heat 
bath (RHB) method to the study of sphaleron transitions in AHM. 

Earlier real-time simulations of AHM ||, [| found that at low temperatures (about 
0.1 of the sphaleron energy) the temperature dependence of the rate qualitatively agrees 
with that predicted by the 1-loop (TST) calculation of Ref . || . The temperature range 
of our simulation is wider and includes somewhat higher temperatures, up to about 1/3 
the sphaleron energy, for which the TST result may no longer be reliable. This allows 
us to estimate the temperature at which TST loses validity. As will be demonstrated 
in the following, this is the temperature at which the scalar field correlation length falls 
below the linear size of the sphaleron. It is also interesting to determine the sign and 
magnitude of the rate deviation from the TST prediction. If the rate we measure falls 
considerably short of the latter, it might be indicative of the entropic rate suppression 
which reflects the difficulty of creating a coherent configuration in high-temperature 
plasma. Our results, presented in the following, do indeed show dramatic slowdown of 
the rate growth. 

Our starting point is the (1+1) -dimensional lattice AHM Lagrangian which in suit- 
ably chosen units reads M 



M 4° - 4° 



J 



a 2 |0i+i - exp (iaA] d+1 ) <^| 2 - - (|<^-| 2 - l) 



(1) 



Here j labels sites of a chain whose lattice spacing is a. The temporal component of a 
vector potential, A®, and a complex scalar field <pj reside on sites of the chain, whereas 
the spatial component of the vector potential, resides on links. Imposing the 

A = condition one obtains a Hamiltonian (we shall drop the Lorentz index of A 1 from 
now on) 



2 z ^ 
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where ?Tj and -Ejj+i are canonically conjugate momenta of <pj and Ajj +1 , respectively. 
The Hamiltonian equations of motion 



3 



E jd+1 = i4>j (exp(iaA jd+1 )(j)* +1 - 0*J + h.c, 
1 . 



a J 



■kj = - (exp(iaA jtj+1 )<f)* j+1 + exp(-iaA i _ 1>J )0*_ 1 - 



a 



-a0*(|0/-l) (3) 
are supplemented by the Gauss' law constraint 

~ (Ej, j+ i - Ej-v) = Im (tt^*) . (4) 

The same dynamics can be described in terms of gauge-invariant variables. To this end, 
the Higgs field is rewritten in polar coordinates: <pj = p.,- exp(i<x,). Defining = 
otj+i — otj — aAjj + i, €j t j+i = ^Ejj + i and introducing canonical momentum %j for pj we 
see that (||) together with (|) is equivalent to 

CI 

i _ 1 ( e j+l,j+2 ~ <.,..j l ' j ri t c 



Pj 



a \ Pj+i Pj 
1 



-7T P 

a J 



7r| = 3 ° — + ^ (Pj+l COS ftjj+i + Pj_i COS fej-ij - 2pj) 



a p 3 j 

a Pj (pf - l) . (5) 



The equations of motion in this form involve only two pairs of real canonical variables, 
namely, pj, and ejj+i, fyy+i- It is easy to see that fl5|) follow from the Hamiltonian 

i 

obtained from polar-coordinate form of (|2]) by substituting (|4]). 

In the following we shall use both presented forms of the equations of motion. On 
one hand, the Cartesian form (0) allows better numerical handling of sphaleronlike field 
configurations in which the Higgs field is close to zero at one or more sites. For this 
reason we use it for real-time evolution in the bulk of an open gauge-Higgs system. On 
the other hand, the gauge-invariant form (|5|) involves less degrees of freedom and lends 
itself easier to linearization. We therefore use it for the heat bath construction. 

As a heat bath we take AHM linearized in the vicinity of one of its gauge-equivalent 
vacua, which, as is well known, is a system of two free fields: the radial Higgs field g 
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and the gauge field e whose masses are a/2 and respectively. Those are coupled at 
the boundary site of the AHM ((||) or (|)) each to its interacting counterpart, i.e. g to 
p, and £ to e. Suppose for definiteness that the left boundary of the interacting system 
separating it from the linear heat bath is at j = site of the chain. The field equations 
at the boundary are then modified compared to (|^), namely, the second and the fourth 
equation of @ are replaced by 

6-1,0 = i(eo,i-2*e_ li o)-a^e_ li o + J D( V / e,[e-i,o])+i 71 ( v // e,t), 

n p = -(p 1 cos6 ,i-2po + l)-2a(p -l) + J D(v / 2,[po-l]) + i 71 (v / 2,t). 
d 

(7) 

In going from ([|) to (0) we linearized in the vicinity of po = 1, 6_x,o = and, following 
Ref. ||, introduced two terms describing interaction with the linear heat bath at the 
boundary. The D(m, [a]) term represents the reaction to the motion of a boundary field 
a from the heat bath. The mass of the corresponding heat-bath field is m. Explicitly, 

D(m, [a]) = f a(t') Xm {t-t')dt\ (8) 

J — CO 

where the Fourier image of the causal response function Xm(t) is 

X m (uj) = 2isign(a;)y / (a; 2 - m 2 )(l + a 2 (m 2 - uj 2 )/A) (9) 

for frequencies m < \u\ < ^m 2 + A/a? and vanishes outside this range. The D(m, [a]) 
term thus describes dissipation processes: the choice of \m 

it) ensures that the waves 

traveling across the boundary into the heat bath are completely absorbed. In order for 
the system to reach thermal equilibrium with the heat bath, the F(m, t) term describing 
thermal fluctuations of a heat-bath field at the boundary should be included. According 
to fluctuation-disspation theorem, F(m, t) is a Gaussian random variable whose time 
autocorrelation is related to Xm{t)'- 

(F{m, t)F{m, t + r)) = 9 f Xm(t')dt', (10) 

Jo 

where 9 is the temperature. We solve numerically the system of equations (0) together 
with (0) and its right-boundary analog. Numerical implementation of the boundary heat 
bath is explained in detail in Ref. As has already been mentioned, the Cartesian 
form of equations of motion in the bulk is used due to its superior numerical properties. 
On the other hand, we use polar coordinates to evolve boundary fields. In order to be 
able to transform from polar to Cartesian coordinates at the boundary, we need to keep 
track of the angular variable a of the boundary Higgs field. This is done with the help 
of Gauss' law which for the linearized system gives 

«i = ~ ( E o,j+i ~ E o-^o) ■ (U) 
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We use second-order Runge-Kutta algorithm for numerical solution of (0^,0). A num- 
ber of criteria were applied in evaluating the algorithm performance and in determining 
the value of the time step. In particular, we tested the absorption properties of the simu- 
lated heat bath by dropping the noise term from the equations of motion and cooling an 
initially hot system. The resulting cooling curve was then compared to that of the same 
system in a large real zero-temperature heat bath. In a similar way we studied thermal- 
ization of an initially cold system in the heat bath. The temperature was measured by 
averaging the kinetic energy of the radial Higgs field over the system and over the time 
history. In all our simulations the temperature of a thermalized system was found to be 
within 3% of the assigned value. 

One special numerical issue to be dealt with in a real-time simulation of a gauge 
theory is the accuracy of the Gauss' constraint. The equations we solve are consistent 
with the Gauss' law. However, numerical errors give rise to a small spurious static charge 
density which should be kept in check in order to have no impact on the quantities we 
measure, in particular, the sphaleron transition rate. Since the rate is exponentially 
sensitive to the sphaleron energy, we computed the perturbative correction to the latter 
in presence of a small static charge distribution q(x), AE sph ([q}). The corresponding 
expression is derived in the Appendix. We then used pAE sp h([q]) averaged over the 
sphaleron positions as a criterion of Gauss' law violation (here and in the following j3 
denotes inverse temperature). In all our simulations the amount of the spurious charge 
was far too small to have any measurable impact on the sphaleron transition rate. 

Our attention in this work is focused on the temperature dependence of the rate. 
For this reason we performed all our simulations at a fixed value of the gauge coupling 
£ = 10, with the exception of preliminary study at £ = 0.5. Most of our measurements 
were done for a chain of length L = 100 and lattice spacing a = 0.5. Since our results are 
to be compared to the analytical prediction for AHM in the continuum, we checked their 
dependence on both lattice cutoffs by performing additional simulations at a = 0.25 and 
at L — 50. The Runge-Kutta time step was 0.01 for a = 0.5 and 0.004 for a = 0.25. 
For each set of a, L, and the temperature the simulation time was 10 5 time units. To 
add confidence to our measurements we also performed microcanonical simulations at a 
number of parameter values. The agreement with the canonical results was good except 
for the (3 — 5, L — 100, a = 0.5 case for which the canonical simulation gives somewhat 
higher value of the rate. 

Following Ref. [|J, we extracted the sphaleron transition rate from Acs(i), the time- 
averaged squared deviation of the Chern-Simons variable Nqs = (2vr) -1 / A(x)dx for a 
lag t. For lags shorter than the average time between consecutive sphaleron transitions 
Acs(t) is determined by fluctuations of Nqs in the vicinity of one of its vacuum values. At 
lags much longer than the lifetime of a vacuum many uncorrelated sphaleron transitions 
would have occurred, each changing the value of Nqs by an integer, and a random-walk 
behavior sets in: 

A cs (t)=TLt, (12) 

where T is the sphaleron transition rate per unit length. Figure [I] illustrates the described 
lag dependence of A cs (t). 
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Figure 1: Lag dependence of the average squared deviation of the Chern-Simons variable. 

Once Acs(^) is known, it can be (at large enough values of t) fitted to a straight 
line through the origin to yield T. Since the values of Acs(^) at different t are strongly 
correlated, we found that, while the quality of fit remains high as more and more Acs(^) 
data points are included, the error on T is not reduced significantly by fitting with more 
degrees of freedom. We therefore simplified the procedure and extracted T from a single 
value of A cs (t) at t = 1000. This choice of a lag is suitable since, on one hand, it is 
much shorter than our total simulation time (10 5 ), while on the other hand it is at least 
several times longer than the average time between consecutive sphaleron transitions in 
the temperature range considered. We also verified that T remains constant in a wide 
range of lags including the chosen one. 

Summary of all our rate measurements is presented in Table [l], while for Figure @ we 
selected the results that best reflect the important features of T dependence on the inverse 
temperature /3, as well as on a. Obviously, no measurable dependence on L is observed. 
The absence of finite-size effects is to be expected of our heat-bath construction. Namely, 
at low temperatures the linearized heat bath closely imitates the infinite extension of the 
nonlinear system beyond the boundaries. At high temperatures the system becomes 
less and less correlated in space and time, and, as a result, the boundary effects lose 
importance. We also observe no dependence of the rate on a, except for the highest 
temperature {(3 = 3) considered. The virtual independence of T of the lattice spacing at 
low temperatures has been found in earlier work [^|, and is confirmed by our results. 
As will be shown shortly, the difference in the rates between the a = 0.5 and a = 0.25 
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Figure 2: Temperature dependence of the transition rate T. The solid curves correspond 
to the TST prediction (13), with the value of £ indicated near each curve. 



cases at (3 = 3 is consistent with other properties of the model at this temperature. Our 
rate measurements are to be compared to the TST prediction || 

u_ ( Q(3E S ^ f T(a + s + l)T(a-s) \^ 

v= 2^\-^r) { r(« + i)r(a) J ex PH^)> ( 13 ) 

where the sphaleron energy _E sp h = 2y2/3, a 2 = s(s + 1) = 2£, and u) 2 _ = s + 1 is the 
negative squared eigenfrequency corresponding to the sphaleron instability. 

As Figure |] clearly shows, the values of T at low temperature are close to those given 
by (p~3|) . While the approach of measured T to the TST prediction is slow, the two nearly 
coincide at the lowest temperature considered, (3 = 14. This nice agreement shows once 
again that our heat bath construction works as intended. But the most notable feature of 
our results is their dramatic departure from the TST-predicted values starting at about 
(3 = 5. For £ = 10 the discrepancy is a factor of 5 already at (3 = 5, and grows at (3 = 3 
to a factor of 10 for a = 0.5 and a factor of 20 for a = 0.25. The situation is similar for 
£ = 0.5. Moreover, in the latter two cases the rate practically does not grow between 
(3 = 5 and (3 = 3. Note that the deviations are much larger than our measurement error 
bars and are therefore statistically significant. 

It is natural to ask why this strong lagging of the measured rate behind the TST one 
begins in the vicinity of (3 = 5. To this end recall the underlying assumption of 
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2.89 


0.09 ±0.01 
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8.87 


(37 ±4) x 10- 
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10 


0.5 


100 


3.91 


0.068 ± 0.007 




10 


0.25 


50 


2.99 


0.046 ± 0.004 




10 


0.5 


100 


4.95 


0.034 ±0.003 




10 


0.25 


50 


4.99 


0.037 ±0.004 




10 


0.5 


100 


6.89 


0.015 ±0.002 




10 


0.25 


50 


7.02 


0.016 ±0.002 




10 


0.5 


100 


8.90 


(37 ±4) x 10" 
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10 


0.25 


50 


8.98 


(33 ±5) x 10" 


-4 


10 


0.5 


100 


11.85 


(32 ± 3) x 10" 
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0.5 


0.5 


100 


2.97 


(78 ±8) x 10" 


-4 


10 


0.5 


100 


13.86 


(60 ± 7) x 10- 
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0.5 


100 
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(67 ±7) x 10" 
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50 
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0.033 ±0.005 
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50 
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0.016 ±0.002 
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100 


8.95 


(54 ±7) x 10- 
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10 


0.5 


100 


4.94 


0.055 ±0.007 




10 


0.25 


50 


3.00 


0.040 ± 0.004 




10 


0.5 


100 


8.95 


(40 ± 5) x 10" 


-4 















Table 1: Summary of transition rate measurements. The inverse temperature (3 is given 
as deduced from the average kinetic energy of the radial Higgs field. 



Chern-Simons number diffusion is dominated by evolution of configurations resembling 
the vacuum into those resembling the zero-temperature sphaleron. It is clear that the 
Ncs diffusion can only be described in these terms as long as the Higgs field is correlated 
on a length scale larger than the sphaleron size (2\/2)- The corresponding correlation 
length A can be found by measuring a gauge-invariant two-point function [[0| 

Cji = 4>*(j)i exp -ia A k ,k+i = PjPi exp -i ^ h,k+i • (14) 
V k=j J V k =i ) 

A rough estimate of A at low temperatures is obtained by averaging Cji over the thermal 
ensemble with H' of @ replaced by its linearized version. Performing Gaussian integra- 
tion over the b variables one finds A = 2(3. This is, in fact, an overestimate of A, since 
thermal fluctuations of the radial Higgs field are not taken into account. Figure |3] shows 
the values of A obtained by fitting (Cji) to const x exp(— \j — /|/A). As expected, the 
breakdown of the saddle-point approximation for the rate occurs as A becomes smaller 
than the sphaleron size. Note that this is true for both values of £ considered. At /3 — 3 
A = 1.47 for £ = 10, only about 3 times larger than the lattice spacing a = 0.5. Hence the 
field strongly fluctuates at length scales comparable to the lattice spacing. It is therefore 
not surprising that we find the a dependence of the rate at this temperature. 

At this point it is unclear what causes the sharp slowdown of the rate growth at 
the high-temperature end of our measurement range. We cannot exclude a possibility 
that at temperatures in question crossing the Acs — half — integer separatrix in the 
configuration space of the model JTJ in close vicinity of the sphaleron saddle point is 
still strongly preferred energetically, but is already suppressed entropically. It is usually 



assumed [|| y, 0, |15| that at temperatures above the sphaleron energy the rate grows 



like a power of the temperature. It could be that what we observe at (3 < 5 is a crossover 
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Figure 3: Temperature dependence of the correlation length deduced from the gauge- 
invariant two-point function (14). The sphaleron size is shown by the dashed line for 
comparison. 



from the exponential to power-law behavior of the rate. It is not clear, however, that 
at such a crossover the rate should stop growing as it does for £ = 10, a = 0.25 and for 
£ = 0.5, a = 0.5. The only way to resolve this puzzling situation is by rate measurements 
at still higher temperatures, as well as smaller lattice spacings. We plan to do so in the 
future. 

To summarize, we performed an accurate measurement of the sphaleron transition 
rate in AHM averaged over the canonical ensemble. The latter was obtained by im- 
mersing the system in a realistic heat bath. The ergodicity of real-time evolution was 
thus achieved without having to introduce an artificial viscosity parameter. Our rate 
measurements approach the corresponding TST estimate at low temperature. This is in 
agreement with Ref. || where similar measurements were performed microcanonically. 
The highest temperature studied in that work was 0.103-Esph, well within the range of 
applicability of TST. In going beyond that range, we found dramatic slowdown of the 
rate growth, with suppression factor as large as 20 relative to TST at the highest tem- 
perature considered. Our measurements show that the breakdown of TST occurs as soon 
as the correlation length deduced from Cji flU])) becomes comparable to the sphaleron 
size. This suggests that a similar object might serve as a criterion for applicability of the 
sphaleron approximation in other theories, including the realistic 3+1- dimensional case. 
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Appendix 

In this Appendix we give a perturbative estimate of the shift in the sphaleron energy 
in presence of a small static charge density q. For simplicity we use continuum, rather 
than lattice, formulation of AHM. We shall also assume that the system has an infinite 
length. Analogous to (0), the Hamiltonian depends on three pairs of canonical vari- 
ables A, E, p, 7c r ho, a, 7T a . Static configurations are obtained by minimizing the energy 
with respect to all the variables on a subspace constrained by Gauss' law 7r a = E' — q 
(prime means derivative with respect to the spatial variable x). Vacuum configurations 
correspond to the absolute minimum of the energy on that subspace, while sphalerons 
minimize the energy among configurations with p = at one point. Two of the vari- 
ables, A and a, enter the Hamiltonian only in combination b = a' — A, thereby effectively 
reducing the number of variables by one. It is convenient to define a new complex field 

$(x) = p(x) exp (^i J b(x')dx'^j . (15) 

After eliminating n a with the help of Gauss' law one obtains 



H = ~ 1 <lv 



(16) 



£E 2 + vr 2 + {E 'f + |$'| 2 + ^(|$| 2 - l) 2 



Extremization of energy leads, apart from the trivial condition tt p = 0, to a system of 
coupled equations for <E> and E: 



E'-q 

|$|2 



0,$" + -^— ^-$(|$| 2 -1) = 0. (17) 



These equations are simplified if we introduce an electrostatic potential for the electric 
field: E = Y'. If we require that q vanishes for \x\ above certain value, E and |$| must 
respectively approach and 1 as x — > ±oo. The electrostatic potential Y will then 
approach a constant value. If that constant is chosen to be 0, (fi~7|) takes form 



Y" - £|$| 2 r = q; $" + £ 2 F 2 $ - $(|$| 2 - 1) = 0. (18) 

We are interested in perturbative corrections to the lowest order in q to the vacuum 
Y = 0, $ = 1 and to the sphaleron configuration Y = 0, $ = tanh(x/v / 2)- We will 
show in the following that in both cases the correction to Y is first order in q. It is then 
clear from the second equation of flI%D that the correction to $ is at best second order. 
Inspecting the Hamiltonian (|T6|) and bearing in mind that the unperturbed solutions 
extremize H for q = we conclude that the correction to the energy is second order in 
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q and comes solely from the first-order correction to Y. Our task therefore reduces to 
solving the first equation of fll8D in the unperturbed $ background. Writing Y(x) as 
/ G(x,y)q(y)dy we obtain an equation for the Green's function G(x,y): 

(d 2 x -^)\ 2 )G(x,y) = 6(x-y). (19) 

It is easy to verify that 

G v (x,y) = -^=exp (-y/t\x -y\j (20) 

for the vacuum (<£>(x) = 1). For the sphaleron configuration (3>(sc) = tanh(x/v^2)) the 
solution of (^) is also straightforward but somewhat cumbersome. The result can be 
expressed in terms of hypergeometric functions: 



G s (x,y) = -Ccosh-v^x/v^coshrv 7 ^ 

exp (x/ \pl 



xF «+,«_; 7; 



2 cosh (x/y/2) 
exp (-y/y/2) \ 

* Fia+ ' a --' r ' 2<l(y,4 ) (2i; 



if y > x, with x and y interchanged otherwise. Here 



1 + ± y/l+M l + a+ + a. 



a ± = ; 7 

C = r ^ +)r(a " } . (22) 

4 1+ v 2 «r( 7 )r( v / 2e) 

Substituting the correction to electric field into (|1^) we find the energy shift of a state 
due to the static charge q: 

AH = -^J dxdyG(x, y)q(x)q(y). (23) 

This result is hardly surprising: G(x,y) is nothing but the Coulomb potential at x due 
to a unit charge at y. Therefore to the lowest order in q the energy shift is equal to the 
Coulomb energy of external charge distribution. The correction to the sphaleron energy 
barrier then follows immediately: 

A£ sph = -- J dxdy (G s (x, y) - G v (x, y)) q(x)q(y). (24) 

References 

[1] M. Shaposhnikov, Nucl. Phys. B (Proc. Suppl.) 26 78 (1992). 



12 



[2] For a review see, e.g., A. D. Dolgov, Phys. Rep. 222, 309 (1992); M. Shaposhnikov, 
Physica Scripta T 36 (1991). 

[3] J. Ambjorn, T. Askgaard, H. Porter, and M. S. Shaposhnikov, Phys. Lett. B 244 
479 (1990); Nucl. Phys. B 353 346 (1990). 

[4] A. I. Bochkarev and Ph. de Forcrand, Phys. Rev. D 44 519 (1991). 

[5] A. Krasnitz and R. Potting, A realistic heat bath: theory and application to kink- 
antikink dynamics, IPS Research Report 92-30 (1992), to be published in Nuclear 
Physics B; Nucl. Phys. B (Proc. Suppl.) 30 900 (1993). 

[6] A. I. Bochkarev and G. G. Tsitsishvili, Phys. Rev. D 40 1378 (1989). 

[7] M. Alford, H. Feldman, and M. Gleiser, Phys. Rev. Lett. 68, 1645 (1992). 

[8] A. I. Bochkarev and Ph. de Forcrand, Nonperturbative evaluation of the diffusion 
rate in field theory at high temperatures, IPS Research Report 92-18 (1992); Nucl. 
Phys. B (Proc. Suppl.) 30 719 (1993). 

[9] D. Yu. Grigoriev, V. A. Rubakov, and M. E. Shaposhnikov, Nucl. Phys. B 326 737 
(1989). 

[10] D. Yu. Grigoriev and M. E. Shaposhnikov, Int. J. Mod. Phys. A 5 4677 (1990). 

[11] D. Yu. Grigoriev, M. E. Shaposhnikov, and N. Turok, Phys. Lett. B 275 395 (1992). 

[12] W. H. Tang, Simulation of the topological transitions for the (1+1) -dimensional 
Abelian Higgs model at high temperature, University of Amsterdam thesis (1992). 

[13] E. Mottola and A. Wipf, Phys. Rev. D 39 588 (1989). 

[14] J. Kripfganz and C. Michael, Lattice study of sphaleron transitions in a 2D 0(3) 
sigma model, University of Liverpool research report LTH 298 (1993). 

[15] M. Dine, O. Lechtenfeld, B. Sakita, W. Fischler, and J. Polchinski, Nucl. Phys. B 
342 381 (1990). 

[16] Da-Shin Lee and D. Boyanovsky, Dynamics of phase transitions induced by a heat 
bath, University of Pittsburgh research report PITT-93-4 (1993). 

[17] J.S. Langer, Annals of Physics 54 258 (1969). 

[18] For a recent review on transition-state theory see, e.g., P. Hanggi, P. Talkner, and 
M. Borkovec, Rev. Mod. Phys. 62 251 (1990). 

[19] J. Bricmont and J. Frohlich, Phys. Lett. B 122 73 (1983). 



13 



